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Complex asteroseismology of the hybrid B-type pulsator 7 Pegasi: 
a test of stellar opacities 
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Using the updated oscillation spectrum of 7 Pegasi, we construct a set of seismic models which reproduce two pulsational 
frequencies corresponding to the ^ = 0, pi and I = 1, gi modes. Then, we single out models which reproduce other 
well identified modes. Finally, we extend our seismic modelling by a requirement of fitting also values of the complex, 
nonadiabatic parameter / associated to each mode frequency. Such complex asteroseismology of the B-type pulsators 
provides a unique test of stellar metallicity and opacities. In contrast to our previous studies, results for 7 Peg indicate that 
both opacity tables, OPAL and OP, are equally preferred. 
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1 Introduction 

Oscillation spectra of the main sequence B-type pulsators 
are characterized by a low number of frequency peaks and 
a lack of distinct regularities. Despite of these facts, aster- 
oseismology of these massive stars provides valuable con- 
straints on stellar physics and evolution. 

Last years, the most attractive targets for asteroseismic 
studies have become hybrid pulsators which exhibit low or- 
der p/g modes as well as high-order g-modes. This is be- 
cause these two types of modes sound various parts of the 
stellar interior. In the case of the early B-type stars, modes 
typical for the (3 Cep and SPB type occur simultaneously. 
The analysis of extensive observations, especially from mul- 
tisite campaigns, allowed to discover several pulsating vari- 
ables of the P Cep/SPB type. 

The first detecte d hybrid pulsator of the early B spectral 



ine nrst aetecte d nybria pulsator 01 tne early a spectral 
type was v Eridani jAerts et al.ll2004t Handler et al.ll2004: 



Jerzvkiewicz et al. 20 05b. The next exam ple is 12 Lacertae 
(iDesmet et al., ,2009; iHandler et al. I l2006h . Many excellent 
papers devoted to seismic modeling of these stars have been 
published, but a limited space of this article does not allow 
to mention them. 

7 Pegasi is one more hybrid pulsator of early B spectral 
type. Recent analysi s of the space based observations from 
the MOST satellite handler et al. "2009*) and ground based 
photometry and spectroscopy ([Handler 2009) led to discov- 
ery of 14 pulsational frequencies with 8 of the (3 Cep type 
and 6 of the SPB type. 

Asteroseismic modeling, consisting in fitting theoretical 
and observational values of pulsational frequencies, can be 
extended by adding another seismic tool associated to each 
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mode frequency. Such a tool is the amplitude of the bolo- 
metric flux perturbation to the radial displacement, called 
the /-parameter. This seismic probe was introduced by 



Daszvi iska-Daszkiewicz. D ziembow ski & Pamvatnvkhl (12003 
2005). Theoretical values of / are obtained from stellar pul- 
sation computations and their empirical counterparts are de- 
rived from multicolour photometry. The parallel fitting of 
pulsational frequencies and corresponding values of the J- 
parameter was termed complex asteroseimology by 
Daszviiska-D aszkiew icz & Walczak ( 2009 ^. The hybrid pul- 



sators are of particular interest because the dependence of / 
on pulsational frequency, v, and the mode degree, I, is com- 
pletely different for p modes and high-order g-modes. 

Complex seismic modeling has been applied to the 
Cep star 6 Oph (Daszvnska-Daszkiewicz & Walc zak 2009 ) 
and to V Eri (IDaszyriska-Daszkiewicz & WalczakI 2010l) . In 
the case of 9 Oph, we got a strong preference for the OPAL 
data. From the v Eri analysis a contradictory result was ob- 
tained: the /3 Cep-type modes indicate the OPAL opacities 
whereas the SPB-type modes prefer the OP data. 

In this paper we present complex seismic modeling of 
7 Peg. In Section 2, we give a short description of the star. 
Section 3 is devoted to mode identification for all detected 
pulsational frequencies. Results of our seismic modeling are 
presented in Sections 4. Conclusions end the paper. 



2 7 Pegasi 

7 Peg (HD 886) is a bright star (V=2.83 mag) of the B2IV 
spectral type. Its low amplitude variability was discovered 
almost a cen tury ago from the radial velocity observations 
Burnsll911 ). These changes were confirmed by McNamara 



1953b who classified 7 Peg as the /3 Cephei pulsating vari- 



able. 
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For dozens of years, 7 Peg has been thought to be a 
monoperiodic star, pulsating in the radial mode w ith the pe- 
riod of about 0.1517 5 d ( Jerzykiewicz 197d) . Finally, 
Chapellier et al. ( 2006h reported discovery of three additional 
pulsational frequencies; on e in the Cep range an d two in 
the SPB frequency domain. Chapellier et al. ( 20061) claimed 
also that this star is a spectroscopic binary with the orbital 
period of 370.5 d. However, on the basis of space obser- 
vations from the MOST mission and the ground ba sed pho- 



tometry and spectroscopv. ' Handler et al. ( 2009h and Handler 
(12009) showed that 7 Peg is a single star and the hypothet- 
ical orbital variations can be explained by the high-order 
g-mode pulsation. 

In Fig.[T] we show a position of 7 Peg in the HR 
diagram. We included the whole range of the effective 
temp erature, Tpff , available in the Uterature 



2009 



Daszvriska-Daszkiewicz & Niemczura 2005: Handler 



Huang & Gies 20081: iMorel et all I2OO6I) and the lu 



minosity was calculated from the Hippar cos paral l ax an d 
the bolometric correction from models of iKurucj (120041) . 
The evolutionary tracks were computed with the Warsaw- 
New Jersey ev olutionary code adopting the OPAL opacities 
(llglesia s & Rogers J 996.) and the newest heavy element mix- 
ture bv lAsplund et al.l ilO^ . hereafter AGSS09. Only evo- 
lution on main sequence is shown. The star is a slow rotator 
with Vi-ot sin i ~ as determined from the Silll lines by 
Telting et al. (l2006l) . The rotational splitting from the possi- 
ble two c omponents of the £ = 1 triplet gave Vi-ot ~ 3 km/s 
([Handler et al..,2009.) . Metallicity o f 7 Peg is smaller than 



solar value. Morel et al. (l2006h determined 



the 

Z = 0.009 ± 0.002 from the optic al spect ra, whereas 
Daszviiska-Daszkiewicz & Niemczural (l2005b obtained 
[m/H] = -0.04 ± 0.08 (equivalent to Z = 0.018 ± 0.003) 
from the lUE ultraviolet spectra. The evolutionary tracks in 
Fig. 1 were computed at the assumption of the equatorial 
rotational velocity of 3 km/s, hydrogen abundance of X = 
0.7, two values of metallicity parameter Z = 0.010 and 
0.015 and no overshooting from a convective core, aov = 
0.0. Lines labeled as n = 1 and n = 2 will be discussed 
later on. 



3 Identification of the mode degree, £ 

To determine the mode degree, £, of the observed frequen- 
cies of 7 Peg, we made use of the light variations in the 



locity changes ([Handler et al .120091) . To this aim we applied 



Stromgren ttf w passbands (|HandleJ2 009) and the radial ve- 

39|).' 
three approaches. 

In the first case, we compared theoretical and observa- 
tional values of the amplitude ratios and phase differences 
between the available passbands and relied on the theoreti- 
cal values of the /-parameter. The second approach includes 
also the radial velocity variation. In the third method we 
used amplitudes and phases themselves and the /-parameter 
was determined from the observations together with the 
mode degree, £ 
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Fig. 1 The observational error box of 7 Peg in the HR di- 
agram. The evolutionary tracks were computed for two val- 
ues of metallicity, Z. Lines of the constant period (0.15175 
d) for the fundamental (n — 1) and first overtone (n = 2) 
radial mode are also drawn. 



( Daszvriska-Daszkiewicz. Dziembowski & Pamvatnvkhl2003 



20051) . In the case of B-type pulsators, the last method de- 
mands the radial velocity measurements to get a unique iden- 
tification of £. 

In Table[T] we give the most probable values of the mode 
degrees for the 7 Peg frequencies from these three approa- 
ches. The dominant mode, j^i, is certainly radial. The iden- 
tification of 1^2 is also unique; it is undoubtedly the £ ~ 2 
mode. In the case of we did not get a consistent result; it 
can be £ = 1 or ^ = 2. The 1/4, vq, v% and 1^12 modes have 
most probably the degree £ — 2. The vc, and vn frequencies 
seem to be the dipole modes. Identification of the remaining 
frequencies are ambiguous, although in the case of vu only 
higher degrees are possible, with £ = 6 as most probable. 
It is interesting that vu is close to the first overtone radial 
mode in the 7 Peg models, but this identification is excluded 
by our analysis. 

In Fig.[T] we plot lines of a constant period (0.15175 d) 
corresponding to vi for the fundamental and first overtone 
radial mode. The discrimination of the radial order, n, can 
be done from a comparison of the empirical and theoreti- 
cal values of the /-parameter. This is shown in Fig.|2l In 
the top panel, we plot the empirical and theoretical values 
of the /-parameter on the complex plane assuming that the 
dominant mode is fundamental. The theoretical values of / 
were calculated for three different metallicities, Z = 0.010, 
0.012, and 0.015, overshooting parameter aov = 0.0, hy- 
drogen abundance X = 0.7 and the OPAL data. In the bot- 
tom panel, we plot the same but we assumed that vi is the 
first overtone mode. In both cases, we included models that 
are inside the observational error box of 7 Peg. As we can 
see, an agreement between empirical and theoretical values 
of the /-parameter can be achieved only if vi is the radial 
fundamental mode, n — 1, and Z is in the range (0.012, 
0.015). 
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Fig. 3 The overshooting parameter, ctov, as a function of metalHcity, Z, for the seismic models of 7 Peg found from 
the fitting of the V\ frequency (the £ = 0, pi mode) and the V'^ frequency (the ^ = \, g\ mode) for hydrogen abundance 
X — 0.7, the OPAL (left panel) and OP opacities (right panel). Hatched areas indicates models fitting /-parameter for v\ 
and v^. 



Table 1 The most probable identification of ^ for the pul- 
sational frequencies of 7 Peg from three approaches. 



frequency 
[c/d] 



phot, 
theoretical / 



phot.+Kad 
theoretical / 



phot.+Vrad 
empirical / 



=6.58974 


^=0 


£=0 


£=0 


i/2=0.63551 




£=2 


£=2 


i/3=0.68241 


£=l 


1=2 


£=1,2 


i/4=0.73940 




£=2 


£=1,2 


i/5 =6.01616 


£=l,3 


1=1,2 


£=1,0,3 


i/6=0.88550 


£=2 


t=2 


£=2,5 


i/7=6.9776 


1=1 


fcl,2,3,6 


£=? 


i/8=0.91442 


£=2,3,4 


£=\,2 


£=2,3,5 


!/9=6.5150 


£>l 


£=\,2,3,6 


£=? 


i/io=8.1861 


1=1 


£=1,2,3,6 


£=? 


1^11=0.8352 


£=\,2 


£=\ 


£=1,2 


i/12 =6.0273 


1=1 


1=1,2 


£=2,5 


!/i3=9.1092 


t=l 


£=1,2,6 


£=? 


i^i4=8.552 


£>4 


£=6 


£=? 



4 Complex asteroseismology 

To compute p ulsational models, we used the linear nonadi- 
abatic code of Dziembowskil ( 1977 ). We started our seismic 
modeling by fitting two p-mode frequencies: I'l ~ 6.58974 
c/d and = 6.01616 c/d. A survey of pulsational mod- 
els showed that if i/i is the radial fundamental mode, then 
j^5, identified as £ = 1, can be only gi. Here, we assumed 
that the azimuthal number of 1/5 is m = 0. In all compu- 
tations, we assumed the rotation velocity of Kot = 3 km/s 
and adopted the latest determination of the chemical com- 
position AGSS09. In the left panel of Fig. 3, we present our 
seismic models fitting and on the aov vs. Z plane for 
the hydrogen abundance of X — 0.7 and the OPAL opaci- 
ties. We depicted the instability borders with the thick solid 
line for the radial mode, vi, and with the thick dashed line 



for the dipole mode, 1/5. Models that are located below the 
77 = lines are unstable. We plotted also the lines of con- 
stant m asses and eff ective temperatures. Results for the OP 
tables ( Seatonll2005 ) are quite similar, as can be seen in the 
right panel of Fig. 3. In both panels, we marked also models 
which reproduce, within the observational errors, the em- 
pirical values of / (hatched areas) of the radial fundamental 
mode (labeled as I'l) and of the dipole mode gi (labeled 
as 1^5). The region covering models fitting the /-parameter 
corresponding to 1/5 is quite large because of larger observa- 
tional errors in amplitudes and phases. Unfortunately, there 
is no model fitting the /-parameter for i'l and 1^5 simultane- 
ously, regardless of which opacity tables are used. Chang- 
ing hydrogen abundance did not help in achieving an agree- 
ment. 

Let us now consider other well identified modes: 1/2, va, 
vq, vii and vu. If we assume that these modes are axisym- 
metric (m = 0), then we can search for pulsational models 
reproducing additionally one of these five frequencies. In 
the case of high-order g-modes (v2, va, vq, vh), we had to 
include more than one radial order In the considered range 
of the metallicity and overshooting parameter, the V2 fre- 
quency is £ = 2, g22, g23 or g24 mode. The frequency is 
the £ = 2, gig or g2o mode and vq can be the £ = 2, gig 
or gi7 mode. We depicted also the £ = 1, gg mode corre- 
sponding to vii, and the £ = 6, g2 mode corresponding to 
1^14. Results are presented in Fig.|4] As we can see, models 
reproducing vn cross the hatched areas described above. 
This means that there are models fitting three frequencies: 
vi, vc,, vii and the /-parameter for the radial fundamental 
or dipole gi mode. Lines of the V2 and modes run through 
the hatched region of 1/5 . 

Although we managed to derive the empirical values of 
the /-parameter for the high-order g-modes, we could not 
find models reproducing simultaneously the real and imag- 
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Fig. 2 Comparison of empirical and theoretical values of 
/ for the dominant frequency, i^i, on the complex plane. The 
top and bottom panels correspond to the hypothesis of the 
fundamental and first overtone radial mode, respectively. 



inary part of / for any of the SPB-type mode, neither with 
the OPAL nor OP opacities. 

5 Conclusions 

The aim of this paper was to construct seismic models of 7 
Peg which reproduce both pulsational frequencies and cor- 
responding values of the nonadiabatic complex parameter 
/. Although we did not fully succeed, we have showed di- 
rections and problems that need to be solved. 

The first problem we encountered, was that there was no 
seismic model reproducing simultaneously the /-parameter 
of the radial (i^i) and dipole (1^5) modes. Secondly, there was 
no seismic model fitting empirical values of the /-parameter 
for any high-order g-mode. 

These inconsistencies can be caused either by the under- 
estimated errors or/and indicate that some additional effects 
should be included in pulsation modeling. One obvious so- 
lution could be inadequacies in the opacity data. Recently, 
Zdravkov & Pamvatnvkh (l2009h tiave suggested increasing 
opacity by 20-50% around the Z-bump and DOB (Deep 
Opacity Bump) to explain the observed frequency range of 
7 Peg. Our further studies will show whether this can help 
also to solve problems reported in this paper. 




0.006 0.010 0.014 2 0.018 0.022 0.026 

Fig. 4 The same as in the left panel of Fig. [3] but we 
marked seismic models fitting also the high-order g-mode 
frequencies 1^2, 1^4, vq, vh and the high degree mode vu. 
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